Muon spin rotation measurements of the superfluid density in fresh and aged 

superconducting PuCoGa 5 
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We have measured the temperature dependence and magnitude of the superfluid density p s {T) 
via the magnetic field penetration depth A(T) in PuCoGas (nominal critical temperature T c o = 
18.5 K) using the muon spin rotation technique in order to investigate the symmetry of the order 
parameter, and to study the effects of aging on the superconducting properties of a radioactive 
material. The same single crystals were measured after 25 days (T c — 18.25 K) and 400 days (T c — 
15.0 K) of aging at room temperature. Penetration depths A(0) = 265(5) and > 498(10) nm are 
derived for the fresh and aged samples, respectively. The temperature dependence of the superfluid 
density p B (T) / p B (0) = A(0) 2 /A(T) 2 is well described in both materials by a model using d-wave gap 
symmetry. Within the context of a strong-coupling, dirty d-wave model, a zero-temperature gap 
value Ao = 3.0fcsT c o is obtained in the fresh sample for a scattering rate V = 0.0057rfcsT c o, which is 
consistent with Abrikosov-Gor'kov (AG) pair-breaking theory. This Ao should be compared to the 
weak-coupling value for a clean d-wave gap Ao = 2.14fcsT c o. In the aged sample the same model 
yields Ao = 2.4fesT c o for T — 0.0107rfcflT c o. This value of F is much less than required by the AG 
pair-breaking formalism. Furthermore, the aged p s (0) is reduced by at least 70% compared to the 
fresh sample, which is also incompatible with AT c /T c o ~ 20%, according to AG theory. We conclude 
that the data in aged PuCoGas support the postulate that the scattering from radiation-induced 
defects is not in the limit of the AG theory of an order parameter which is spatially averaged over 
impurity sites, but rather in the limit of short-coherence-length superconductivity. We show that a 
model calculation consistent with this assumption fits p s {T) in the aged sample rather well. 

PACS numbers: 74.70.Tx, 74.25.Qt, 76.75.-K 
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I. INTRODUCTION 



Interest in the magnetism and superconductivity of 
/-electron materials has remained strong over the last 
several decades, starting with the discovery- 1 , of heavy 
fermion superconductivity in Ce- and U-based mate- 
rials, and continuing today with the investigation 2 - of 
novel ground states in the tetragonal Ce-based CeTIn5 
(T=Co, tr, Rh) materials. Pu-based compounds exhibit 
a particular richness and complexity as is illustrated 
by the discovery 3 - of superconductivity in PuCoGas at 
an order of magnitude higher transition temperature 
(T c q — 18.5 K) than the structurally similar heavy 
fermion superconductor 2 CcCoIn5 (T c ^2.3 K). One rea- 
son for this complexity in Pu materials is that their / 
electrons sit at the boundary between more fully local- 
ized (as in Am compounds) and more fully itinerant (as 
in U and Np compounds) behavior. 4 

There is evidence that the superconductivity in 
PuCoGas may be unconventional in nature. Recent nu- 
clear magnetic resonance (NMR) and nuclear quadrupole 
resonance (NQR) studies 5 display an absence of a co- 
herence peak, together with a decreasing Knight shift 
and a T 3 -behavior in the spin-lattice relaxation rate 



1/Ti just below T c , suggesting that the superconduct- 
ing order parameter has d-wave symmetry, as in the 
copper-oxide superconductors^ Calculations 7 ^ of the 
electronic structure of PuCoGas show a Fermi surface 
consisting of several cylindrical sheets, which is favorable 
for rf-wave, spin-fluctuation mediated superconductivity. 
It has been suggested^ that PuCoGas belongs to a class of 
spin-fluctuation induced superconductors, lying between 
the high-temperature, copper-oxide superconductors and 
the low-temperature, heavy-fermion superconductors on 
a plot of transition temperature versus spin fluctuation 
temperature^^ 

Additional measurements which test this supposition 
are, therefore, important. This is particularly true be- 
cause the natural radioactivity of Pu ( 239 Pu half-life = 
24, 000 years) creates lattice defects which scatter elec- 
trons and can create impurity bands which partly obscure 
the signature of a pure d-wave superconductor. Such ef- 
fects are evident in the low-temperature NQR spin- lattice 
relaxation rate where the T 3 temperature depen- 

dence gives way to a linear-T behavior below about 0.4T C . 
As we show, the temperature dependence of the penetra- 
tion depth \(T) in PuCoGa 5 is less sensitive to these 
defects, thus providing a clean test of the gap symme- 
try and of the robustness of the order parameter to the 
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presence of pair-breaking defects. The magnitude of A(0) 
is strongly sensitive to defects, however, a fact which we 
are able to explain in terms of the relatively short super- 
conducting coherence length in PuCoGas. Finally, the 
study of radiation effects in superconductors is of prac- 
tical importance as well, because increases in the critical 
current density can be achieved by the pinning of flux 
lines at defects. 

In this paper we present transverse-field (TF) muon 
spin rotation (/xSR) measurements of the in-plane mag- 
netic field penetration depth A in the same single crys- 
tals of PuCoGas after 25 days and 400 days of aging at 
room temperature. A(T) is determined by the spectrum 
of quasiparticle excitations exceeding the superconduct- 
ing energy gap A(T), and is thus a sensitive measure 
of the temperature dependence of the supcrfluid density 
p s (T) oc A(T)~ 2 , and the gap structure. Some of the data 
after 25 days of aging has been reported previously^ 



II. EXPERIMENTAL DETAILS 
A. Sample preparation and experimental setup 

Two crystals of PuCoGas measuring ~ 5 x 6 mm 2 and 
about 1/2 mm thick were grown from excess Ga flvrs^ 
as flat plates with the c-axis normal to the surface. The 
crystals were encapsulated in a polyimide coating about 
70 /im thick to prevent contamination from particle and 
ejecta emission. The encapsulated crystals were then at- 
tached and sealed under a helium atmosphere inside a Ti 
cell for further protection. The cell was cooled using a 
continuous- flow He cold- finger cryostat. 

The experiments were performed at the M20 surface 
muon channel at TRIUMF, Vancouver, Canada. Muons 
entered the cell through a 50 (im Ti-foil window, with 
their polarization rotated vertically 90°, perpendicular to 
the beam momentum. The applied field Hq was along the 
incoming beam direction in the TF mode, and parallel to 
the crystalline c-axis. A negligible fraction of the beam 
stopped in the Ti window and polyimide coating. 

The initial set of measurements (after 25 days) were 
performed with the sample mounted directly on the Ti- 
cell backing with Hq = 60 mT. The second set of mea- 
surements (after 400 days, the exact time determined by 
TPJUMF's beam-time schedule) was carried out in the 
same geometrical and spin configurations as before, but 
with the sample surrounded by a Ag backing plate and in 
applied fields of 60 mT and 0.3 T. All data with T < T c 
were taken in a field-cooled mode. 



B. Data analysis 

In an ideal experiment the muons stop randomly on 
the scale of the flux line lattice (FLL) spacing. Therefore, 
the muon spin precession signal P(t) provides a random 



sampling of the internal field distribution n(B), 12 

/oc 
n(B)exp(i lf ,Bt)dB, (1) 
-oo 

where n(B) = (S(B — B(r))) r is the spatial average of 
a Dirac delta function, 7^ (= 2vrx 135.53 MHz/T) is the 
muon gyromagnetic ratio and B{r) is the internal field. 
This equation indicates that the real amplitude of the 
Fourier transformed muon precession signal corresponds 
to the spectral density n(B). In general, an asymmetric 
n(B), the so-called 'Redfield pattern', is expected when 
the FLL consists of straight, rigid flux lines and the ratio 
k = A/£ is not too large^ In such a case, n(B) possesses 
a high-field tail which characterizes the flux-line core ra- 
dius (or coherence length £); the second moment ((AB) 2 ) 
of n(B) determines the London penetration depth. 
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FIG. 1: The fast Fourier transform spectra of the /J.SK signals 
in PuCoGas at 3 K under 60 mT, after subtracting the sample 
holder signal, in both the 400 day-old sample and the 25 day- 
old sample, -ffo was determined from a fit to the normal-state 
data. 

The fast Fourier transform (FFT) provides a reason- 
ably good representation of n(B) in the FLL state. Fig.[T] 
shows FFT spectra for the /iSR signals from the sample 
after 25 and 400 days at T = 3 K. Here the background 
signal from the sample cell has been subtracted from the 
time spectrum before carrying out the FFT. Some age- 
related differences in the spectra are apparent. First, 
while the fresh sample exhibits a slightly asymmetric line 
shape and a small negative field shift, the aged sample 
shows a nearly symmetric shape with almost no negative 
shift. Second, the linewidth has narrowed appreciably 
with aging. 

Attempts to fit the fresh sample data using standard 
model a 14 ' 15 which incorporate a broadened (due to in- 
strumental and nuclear line-width effects) Redfield pat- 
tern for a perfect FLL were not successful. This is due 
to the fact that the lineshape is only moderately asym- 
metric, and thus unique determinations of A and the core 
radius, which depend upon having a high-field tail in the 
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FIG. 2: Muon TF precession signal at T = 4 K, correspond- 
ing to the first term in Eq. The background component 
was characterized at times between t = 3 — 10 fis, where the 
sample signal is fully relaxed, and subtracted to show the in- 
homogeneous relaxation produced by the flux lattice in the 
superconducting state. The solid line is a fit for a Gaussian 
P(t) in Eq. ©. 

field distribution, cannot be obtained. Most of the de- 
creased asymmetry in the lineshape is due to the fact 
that PuCoGas is a high k (> 100) superconductor with 
a relatively large penetration depth (see below), but part 
may also be due to distortions of the FLL from radiation 
induced pinning centers, as we discuss below. 

As described earlier, reasonable fits to the fresh-sample 
data were obtained with simple Gaussian or Lorcntzian 
line shapes, and these two functional forms yielded the 
same temperature dependence for the linewidthsJ^ Nev- 
ertheless, in an attempt to account phenomenologically 
for any small asymmetry, we also performed fits using 
a sum of two Gaussians with different centroid frequen- 
cies and widthsJ^ When these widths were convoluted 
into a single width (as in Ref. 1 1 61 ) the same tempera- 
ture dependence and overall lincwidth was obtained as 
for the single Gaussian fits, but with somewhat larger 
(= 10 — 15%) uncertainties. 

Accordingly, the TF precession spectra for both sam- 
ples were fit to the sum of two terms, corresponding to 
muons stopping in the sample and background materials 
(either Ti or Ag), respectively: 

A G z (t) = Acos(ut + (f))P(t) 

+ A h coa(uj h t + cj))exp(-alt 2 /2), (2) 

where A and A^, to and cub are the partial asymme- 
tries and central frequencies for PuCoGas and the sample 
holder, respectively, Aq is the total positron decay asym- 
metry (j4o = A + Ab), and </> is the initial phased Here 
we used P(t) — exp(— a 2 t 2 /2). Fits to the data yielded 
A/Ab = 1/2. The background signal from Ti in Hq = 60 
mT applied field was well characterized by a Gaussian 
relaxation function, Gb(t) — exp(— a 2 t 2 /2), with w 



0.014 /is 1 . The relaxation rate from the Ag backing 
used in the measurements on the aged sample is negligi- 
ble. 

The quality of the Gaussian fits for the 25 day-old sam- 
ple is illustrated in Fig.[2]for data taken at T = 4 K. Here 
the G z (t) data (Eq. @) from t — 3 — 10 fxs were fit sep- 
arately, and this long-time Ti signal was then subtracted 
from the total spectrum, leaving only the signal from the 
sample in the superconducting state. One sees that the 
Gaussian form for P(t) gives a satisfactory fit. Thus, 
we conclude that single Gaussian fits give an acceptably 
accurate measure of the linewidths in the 25 day-old sam- 
ple. Regarding the 400 day-old data, a single Gaussian 
fit is obviously appropriate, as can be seen in Fig. [TJ 

III. RESULTS 
A. Flux Pinning 
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FIG. 3: Temperature dependence of the real Fourier ampli- 
tude [transform of AoG z (t)] in the 400 day-old sample showing 
flux trapping. In the top figure the sample was field cooled in 
Ho = 99 mT. In the lower figures the field has been reduced to 
89 mT with trapped flux giving rise to a resonance at ~ 13.4 
MHz. 

We also performed a series of measurements in the aged 
sample to observe the effects of flux pinning. The sample 
was cooled to 4 K in Ho — 99 mT field and a spectrum 
was obtained. The field was then reduced by 10 mT with- 
out changing the temperature and another spectrum was 
taken. The sample was then warmed in Hq — 89 mT 
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field, and several spectra were accumulated at increas- 
ing temperatures. The results are shown in Fig. [3] In 
the field-shifted spectrum at 4 K one observes an un- 
shiftcd line with an amplitude of 90-100% of the sample 
signal, together with a line at reduced frequency from 
muons stopping in the sample backing material. This 
clearly demonstrates strong pinning, and is consistent 
with the temperature-independent, field-cooled suscepti- 
bilities found in fresh^ and agedi£ samples of PuCoGa 5 . 
As seen in Fig. [31 the trapped flux gradually leaks out 
of the sample as the temperature is raised. Note, how- 
ever, that for the data in Fig. 2J where the field is held 
constant at the field-cooled value, the sample asymmetry 
A is independent of temperature, as expected. 



B. Temperature dependence of a 
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FIG. 4: (a) Temperature dependence of Gaussian muon spin 
relaxation rate a for Ho = 60 mT in the 25 day-old and 400 
day-old samples. The dotted line is a guide to the eye. (b) 
Temperature dependence of a in the 400 day-old sample for 
Ho = 60 mT and 0.30 T. 

The temperature dependences of a for the fresh and 
aged samples for H = 60 mT are shown in Fig. 0Ja). 
In both samples a increases with decreasing temperature 
below T c due to the formation of the FLL. This increase 
occurs below about 18 K in the fresh sample, in good 
agreement with the measured T c . In the 400 day-old sam- 
ple the data are consistent with a decrease in T c of about 
3 K, in reasonable agreement with the radiation- induced 
reduction of T c (ss 0.24 K/month) reported for PuCoGas 
samples of slightly different isotopic concentrations ^ 

As seen in Fig. HJa), there is no temperature depen- 
dence to the normal state values of a (denoted a n be- 



low). Assuming the muons occupy the same sites^S as 
in CeRhlns we estimate a nuclear dipolar linewidth of 
about 0.21 fis^ 1 for both of the suggested (|, |, |) and 
(0, ^,0) sites, which is close to the average of the mea- 
sured values <7 n ~ 0.2 /is -1 . Furthermore, there is little 
change in er n with aging (Fig. 2Ja)), indicating that ad- 
ditional aging has little effect on the muon site(s). None 
of our conclusions regarding the superconducting prop- 
erties of PuCoGas depend on knowing these exact sites, 
however. 

We also performed measurements with H = 0.3 T in 
the 400 day-old sample to investigate the field depen- 
dence of a. In some superconductors a displays a non- 
monotonic field dependence which has been attributed^ 
to motion of the FLL in low fields where the interac- 
tions between the flux lines is weak due to their large 
separation. Figure [3Jb) shows that both the magnitude 
and temperature dependence of a remain essentially un- 
changed for 60 mT < H < 0.3 T in the 400 day-old 
sample. Note that Ho = 60 mT is at least twice the 
lower critical field H cl = $ ln(A/£)/47rA 2 w 0.035 T re- 
ported previously^ and 4 — 5 times that implied by our 
measurements of A(0) in the fresh sample, as discussed 
below. Here $ = hc/2e = 2.07 • 10~ 15 Tm 2 is the mag- 
netic flux quantum. 



C. Temperature dependence of p s and A 

The /xSR linewidth from the vortex lattice a v is ob- 
tained by subtracting the temperature-averaged normal- 
state linewidth cf n from the total linewidth a in quadra- 
ture: ex 2 = (j 2 — er 2 . The temperature-averaged er n = 
0.204(3) /is" 1 in the fresh sample at Hq = 60 mT, and 
0.227(3) us' 1 and 0.234(3) ^s" 1 for H = 60 mT and 0.3 
T, respectively, in the aged sample. 

The penetration depth X(T) can be deduced from 
a v {T) as follows^^ 

cr^s" 1 ) = 48300(l-/i)[l+3.9(l-/i) 2 ] 1/2 A" 2 (nm), (3) 

where h = H /H c2 {0) and H c2 {0) = 74 T is the upper 
critical field in the fresh^ sample, so h <C 1. (i? C 2 is larger 
in the aged sampled) Eq. ^ is valid for an isotropic 
extreme type-II superconductor (A3> £) with a hexagonal 
FLL. The /iSR rate cr v measures the rms width (AB 2 ) of 
the field distribution from the FLL: er v = 7 M (A_B 2 ) 1 / 2 cx 
p s /m*, where m* is the in-plane effective mass. In the 
following we assume that the temperature dependence 
of a v reflects the temperature dependence of p s , i.e., m* 
is temperature-independent. Figs. [5] and [6] display the 
temperature dependence of p s {T) and A(T) for the fresh 
and aged samples, respectively, normalized to their values 
at T = 0. The choice of normalization is described below. 
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FIG. 5: Temperature dependence of (a) the normalized /iSR 
linewidth <x v (T) cx p s {T) and (b) penetration depth A(T) for 
Ho = 60 mT in PuCoGas after 25 days aging. The solid 
lines are fits using the d-wave model in the strong scattering 
limit described in the text for Ao = 2.8,3.0, and 3.2 ksT c o, 
where T c0 = 18.5 K, T c = 18.25 K, and scattering rate F = 
0.0057rfc s T c0 . ov(0) = 1.52(5) us' 1 and A(0) = 265(5) nm. 



D. Possible effects of flux-line lattice distortion 
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FIG. 6: Temperature dependence of (a) the normalized pSK 
linewidth cr v (T) cx ps(T) and (b) penetration depth A(T) in 
PuCoGas after 400 days aging. The solid squares are for 
Ho = 60 mT and the open squares for Ho = 0.30 T. The solid 
lines are fits using the d-wave model in the strong scattering 
limit described in the text for Ao = 2.3,2.4, and 2.5 fcsT c o, 
where T c o = 18.5 K, T c = 15.0 K and fitted scattering rate 
T = 0.0l7rfc s r c0 . ffy(O) = 0.43(1) ps' 1 and A(0) = 498(10) 
nm. The dotted curve is the best fit to the fresh sample data 
in Fig. [5] The dashed curve is for Ao = 2.4fcsT c o and 
T = 0.0747rfc s T c0 . 



As stated above, the pSR linewidth is directly related 
to the penetration depth assuming a perfect FLL. In real 
materials, however, some distortion of the lattice always 
occurs due to slight misalignments and pinning of the 
flux lines at defect centers. Therefore, before discussing 
our data in terms of possible theoretical models, it is im- 
portant to address these issues in PuCoGas, particularly 
because radiation damage produces pinning centers. 

It is generally agreed that if the FLL is stable, that 
is, does not melt or move with temperature, then the 
temperature dependence of a v (T) accurately reflects the 
temperature dependence of p s {T). The small scatter in 
the data observed in Fig. [3] is fully compatible with sta- 
tistical fluctuations (where Gaussian statistics allow a ~ 
37% probability for fluctuations greater than one stan- 
dard deviation from the mean), and with non-statistical 
errors associated with data fitting, and thus is not evi- 
dence for FLL motion. Furthermore, as discussed above, 
both susceptibility and /iSR measurements are consistent 
with strong pinning in PuCoGas. These facts, and the 
absence of a field dependence for a(T) in the aged sample, 
lead to the conclusion that the temperature dependence 
of <t v (T) provides a good measurement of p s (T)/p s (0). 

The second issue concerns the absolute magnitude of 
a v (T). It is evident from Fig. 2] that there is a strong de- 



pression of cr v in the aged sample (« 70%) despite only 
a modest 15-20 % reduction of T c . Conventional pair- 
breaking theory, originally developed by Abrikosov and 
Gor'kov (AG), 22 would predict about a 40% reduction 
in p s for a reduction in T c of 20%i 23 i 24 It is, therefore, 
important to investigate other possible causes for the re- 
duced linewidth besides impurity scattering. These in- 
clude: 1) significant normal-state inclusions in the aged 
sample which might give rise to linewidth narrowing and 
2) the effects of distortions to the FLL caused by pinning 
at defects. 

We first discuss the small probability of large scale 
normal-state inclusions, noting that a definitive answer to 
this question requires a detailed knowledge of the number 
and morphology of the sample defects following room- 
temperature aging, which is not known with certainty. A 
239p u nuc i eus decays into a 5 MeV alpha particle and an 
86 keV recoiling 235 U atom. Most of the damage cascade 
is caused by the recoiling heavy U atom which has a range 
of w 17 nmii£ Lattice vibrations at room temperature 
usually cause some of this damage to be self-annealed* 2 ^ 
Note, however, that vacancy clusters or dislocation loops 
with radii < £ ~ 2 nm retain their superconducting char- 
acter via the proximity effect 3£ Unfortunately, reliable 
TEM studies of aged PuCoGas, which could visualize 
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the damage, are presently lacking. We note, however, 
that the field-shifting experiments shown in Fig.[3]do not 
show any evidence for flux trapped inside normal-state 
inclusions, because the asymmetry of the unshifted line 
is about equal to that measured in a field-cooled state. 
Given the sensitivity of our measurements, we thus es- 
timate an upper limit for the normal-state fraction in 
our aged sample of < 10 %, which is slightly less than, 
but of the same order as, estimated from recent EXAFS 
experiments. 27 Hence, we find little compelling evidence 
for large-scale, normal-state inclusions in the aged sam- 
ple, and dismiss this effect as a significant cause of the 
observed linewidth narrowing. 

We now discuss possible effects of FLL distortion by 
pinning of vortices. Normally one expects that disor- 
der in the FLL will lead to linewidth broadening. This 
is in fact predicted for the random pinning of stiff, rel- 
atively straight vortex linesji 3 - Experimental evidence 
for such broadening has been reported for many su- 
perconductors, e.g., Y(Nio.8Pto. 2)2620,— CeRu2,— and 
YBa 2 Cu306.95r^- On the other hand, to explain the ob- 
served narrowing of the linewidth in some layered high 
k superconductors, a model in which point vortices ran- 
domly displaced within layers with no alignment or corre- 
lation between the layers was developed^ In some cases, 
this could be topologically similar to the random pinning 
of segments of flexible vortex lines. If so, then this ef- 
fect would be most dominant at relatively low applied 
fields, where the interaction between vortices is weak- 
est. The fact that we observe the same linewidth in the 
aged sample at both Hq = 60 mT and 0.3 T seems to 
argue against this scenario. Furthermore, this model 
for extreme linewidth narrowing was developed for a 
highly anisotropic (layered) superconductor. The effec- 
tive mass anisotropy in PuCoGas can be estimated from 
the values 3 ^ of the slopes of the upper critical fields for 
H || c- and a-axes^ These data yield only a small mass 
anisotropy of about 60-70%, and thus PuCoGas is not a 
candidate for such a model. We, therefore, conclude that 
the most likely cause of the sharply reduced linewidth in 
aged PuCoGas is impurity scattering. Note that since 
any FLL distortion is, therefore, assumed to broaden the 
linewidth, the reduction in superfluid density deduced 
for the aged sample is therefore a lower limit; it could be 
greater. 



IV. MODELING THE DATA 

A. Dirty d-wave model 

In a superconductor whose electrons are paired in an 
L = (s wave) orbital angular momentum state, p s and 
AA = A(T) — A(0) are relatively temperature-independent 
below about T/T c = 0.3, reflecting exponentially- 
activated quasiparticle excitations over a superconduct- 
ing gap which is non-zero over the entire Fermi surfaced 
This is clearly not observed in PuCoGas, as seen in Figs. 



[5] and where the behavior of p s and AA in both the 
fresh and aged samples is approximately linear at low 
temperatures. 

A pairing state with relative angular momentum L = 2 
(d wave) has been found to produce p s oc T at low 
temperatures in the clean high-temperature copper-oxide 
superconductors^ We have, therefore, compared the 
data in Figs. [5] and [5] to model calculations for the su- 
perfluid density p s in a d-wave superconductor. The cal- 
culations follow from the standard response formul a 23 ' 34 
for quasi-two-dimensional quasiparticles with a cylindri- 
cal Fermi surface that relates the superfluid density with 
the penetration depth, p s (T) oc 1/A(T) 2 , 
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where uj p is the Drude plasma frequency, c is the speed 
of light, and e m is the impurity renormalized Matsub- 
ara frequency e m — (2m + l^ksT. The impurity 
self-energy is calculated within the framework of the 
T-matrix approximation^ For the clean case this ex- 
pression simplifies to 1/A(0) = u> p /c. Here A(<f>;T) is 
the phenomenological gap function given by A(cf>;T) — 
A Q tanh{b^jT c /T — 1) cos 2(f), and Ao and b are phe- 
nomenological parameters . 36 ' 37 We used b = 1.65, as 
for the case of weak-coupling d-wave superconductiv- 
ity, though our results are not very sensitive to this 
parameter. The relevant scales are set by the mag- 
nitude Ao at T — and the T-dependence near T c : 
A(T) ~ 6A 0V /1 - T/T c . The model calculations include 
the effects of strong impurity scattering to account for the 
radioactive decay of the Pu atom. Although PuCoGa 5 
is not strictly a two-dimensional (2D) superconductor, 
electronic structure calculations predict several cylindri- 
cal Fermi surfaces.^^ Thus, 2D model calculations for 
the temperature dependence of the in-plane A(T) should 
yield qualitatively reasonable comparisons with the data. 
In what follows, we refer to this model as the 'dirty d- 
wave' model. 

The data for the fresh sample are compared to the 
model calculations in Fig. where cr v (0) and T c were 
used as adjustable parameters. Curves for a range of gap 
parameters Ao are drawn as solid lines to show the sensi- 
tivity of the data to the model parameters. Good overall 
agreement with the data is obtained with this model us- 
ing a scattering rate T — 0.0057rfcsT c o, a value consistent 
with the reduction in T c according to conventional AG 
pair-breaking theory (see below). The measured penetra- 
tion depths and derived parameters are shown in Table 
I. It is worth noting that the penetration depth of our 
fresh sample (A(0) = 265(5) nm) is in good agreement 
with NMR measurements (A(0) « 250 nm)^ The pa- 
rameter n in the relation Ao = nksT C Q is an indirect 
measure of the coupling strength of the pairing. For a 
clean weak-coupling d-wave state n = 2.14. Our best fit 
to the fresh sample data yields n = 3.0(1), indicating an 
enhanced coupling strength. 
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TABLE I: Parameters derived from comparing fiSR rates c v (T) in 25 day-old (fresh) and 400 day-old (aged) PuCoGas to a 
model of dirty d-wave superconductivity with gap Ao, penetration depth A(0) (at T = 0), impurity scattering rate V, and the 
transition temperature of the nominally pure sample T c o = 18.5 K. Note that within AG pair-breaking theory the suppressed 
T c — 15.0 K corresponds to V = 0.0747rfesT c o, which is more than seven times bigger than the V needed to fit the T dependence 
of p a - 



Samples 


Tc (K) 


<7 V (0) (/IS" 1 ) 


A(0) (nm) 


n = A /k B T c0 


V/-Kk B T c0 


Fresh 


18.25(10) 


1.52(5) 


265(5) 


3.0(1) 


0.005 


Aged 


15.0(1) 


0.43(1) 


498(10) 


2.4(1) 


0.010 



The temperature dependence of p s {T), or AA(T), ex- 
hibits an essentially linear behavior at low tempera- 
tures. There are at least two phenomena which might 
be expected to destroy this linearity in a d-wave su- 
perconductor. The first is non-local dynamics in clean 
superconductors ; 39 ! 40 which causes a T 3 temperature de- 
pendence in p s below a temperature ksT* ss (£/A(0))Ao. 
We estimate T* ps (2nm/265nm) • 3 • 18.5 K = 0.4 K. 
Thus, we do not anticipate non-local effects to be signif- 
icant in our measurements. 

Strong impurity scattering can induce a T 2 de- 
pendence in AA(T) at low temperatures in a d- 
wave superconducto n 41 i 42 ' 43 i 44 Conventional AG pair- 
breaking theory, originally developed for magnetic im- 
purities in an s-wave superconductor^ has also been 
found to be applicable to non-magnetic impurities in 
a d-wave superconductor^ The reduction in T c for 
AT C /T C <C 1 in the presence of impurity scattering is 
given by AT C = ttT/4. For our fresh sample AT C = 0.25 
K, yielding T = 0.005irkBT c o, the scattering rate used to 
model the data in Fig. [5] The cross-over temperature 
to T 2 behavior for unitary (strong) scattering has been 
estimated 4 ^ (independently from the mode l 36 i 37 used in 
Figs. [5] and EI to be T cr ps 0.83vTA , where A is the 
maximum gap amplitude. Taking Ao = 3ksT c o for the 
fresh sample yields T cr = 3.3 K. The crossover from T- 
linear to T 2 behavior is expected to be smooth, so that 
AA(T) ps bT 2 /(T + T cr )4i Therefore, a T 2 behavior is 
only anticipated for temperatures significantly below T cr ; 
practically, this works out to be T < T cr /3, as observed 
by /zSB4£ for Zn doping in YBa 2 Cu307_(5. Thus, using 
either our model or this simpler analytic estimate, we do 
not expect to see evidence for strong T 2 behavior as long 
as T > 3 K in fresh PuCoGa5, and we do not. 

The same model has been used to fit the aged sample 
data, shown in Fig. [6l where again a range of gap param- 
eters is displayed. The best fits are obtained for n = 2.4 
and r = 0.0l7rfcsT c o. For comparison, the dotted curve 
in Fig. [5] shows the best fit to the fresh sample data 
(Ao = 3ksT C Q, r = 0.0057rfcsT c o), indicating either that 
the coupling strength has been reduced with aging, or 
that pair-breaking has noticeably reduced the gap value 
by about 20% after 400 days of aging. The dashed curve 
in Fig. [6] is a calculation for a dirty d-wave supercon- 
ductor with A = 2Ak B T c0 and T = 0.0747rfc B T' c o in the 
strong scattering limit. This value of T is required by 
conventional AG pair-breaking theory to reproduce the 



reduction of T c in the aged sample, but it clearly leads 
to poor agreement with the temperature dependence of 
the superfluid density. Furthermore, as we now discuss, 
this model does not account for the strong reduction in 
p s in the aged sample. 



B. Short coherence length model 

In our discussion above, we concluded that the most 
likely cause of the sharply reduced linewidth in aged 
PuCoGas is impurity scattering. This situation is not 
unprecedented. In comparing doped and radiation- 
damaged YBa 2 Cu307_ s (YBCO) superconductors to the 
conventional AG pair-breaking theory, it was found in 
many cases that T c was remarkably robust, even though 
p s was easily suppressed i 46 i 47 ' 48 ' 49 i 50 In Ni-doped and He- 
irradiated YBCO materials, for example, a suppression 
of T c by about 20% was accompanied by a suppression 
of p s by about 70%. This large T c /p s ratio is in con- 
tradiction to conventional 'dirty d-wave' theory, and has 
been addressed theoretically by accounting for the sup- 
pression of the order parameter around the vicinity of 
the defectsj 24 ' 51 ! 52 The conventional treatment assumes 
that the spatial variation of the superconducting order 
parameter (of order £) is large compared to the average 
distance I between the scattering centers, but large com- 
pared to the lattice parameter ao, so that the effective 
order parameter is a spatial average, instead of having a 
large suppression only near the impurities. As expressed 
by Franz et al£^, the effect of a spatially inhomogeneous 
gap is enhanced by a short coherence length relative to 
the lattice parameter, £/ao ps 2 - 5j21 As in YBCO, the 
coherence length in PuCoGas is relatively small (~ 2 
nm in both materials) compared to many superconduc- 
tors where conventional AG pair-breaking theory applies. 
Although we do not know the actual spacing I between 
defects after self-annealing, we can estimate the mean 
distance d between Pu atoms which have decayed ran- 
domly after 400 days as a crude measure of the spacing 
between damage cascades. We find d = 20 nm, which 
is ps 10£. Thus, a spatially inhomogeneous gap model 
with a short coherence length, as considered by Franz 
and coworkers j2£ may explain our observations. 

In this regard, we note that both the suppression of T c 
with aging and the value of £ are about twice as large in 
PuRhGa 5 compared to PuCoGas^^ 4 -^ Very recently, 
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an enhanced suppression of the superfluid density with 
chemical doping has also been reported for the related 
compound CeCoIns-^Sn^.— 

In order to test our hypothesis of short-cohcrcncc- 
length superconductivity for PuCoGa 5 , we calculated the 
superfluid density by solving the Bogoliubov-de Gcnncs 
(BdG) equations of a 2D <i-wave superconductor: 

E(^)0tH(|). m 

Here (M™,ff) T are the eigenfunctions at site i corre- 
sponding to the quasiparticle excitation energy E n . The 
normal-state single-particle lattice Hamiltonian is 

Hij = -t5 i+ s,j + {U imp - p)S i3 , (6) 

where t is the hopping integral between a specific lattice 
site and its four nearest neighbors as denoted by 5 — 
(±1, 0) and (0, ±1), p is the chemical potential, and C/; mp 
is the impurity potential modeling the on-site disorder. 
The self-consistency equation for the gap function of a 
d-wave superconductor on a lattice is given by 

A iJ=i+5 = j Y} U >T +U K*] tanh(E n /2k B T) , (7) 

n 

with V being the pairing strength. Note that the quasi- 
particle energy is measured with respect to the chemical 
potential. We follow an iterative procedure to solve self- 
consistently the BdG equations by exact diagonalization 
on a 20 x 20 lattice, and, for comparison, on a 24 x 24 
lattice. Using a suitable guess for an initial Ay, a new 
one is then calculated, and the process is iterated until 
the desired convergence is achieved. 

The linear response calculation of the superfluid den- 
sity p s for the lattice model follows the approach de- 
scribed by Scalapino et al.— and was applied to 2D d- 
wav o 24 i 58 and s-wav o 59 ' 60 superconductors. A detailed 
technical description of the response calculation will be 
given elsewhere*^ Our calculations of the transition tem- 
perature T c , the lattice-averaged zero-temperature or- 
der parameter Ao = (Ay(0)) and superfluid density 
p s (0) are in good agreement with those of Franz et 
al.— for the same set of parameters as used by them 
(V = l.lM,U imp = lOOt, ^ = -0.36*). Because of finite 
size effects of our lattices, we averaged our calculations 
over 20 random impurity configurations at fixed concen- 
tration, as indicated by the statistical error bars in Fig. [7] 

In Fig. [7] we compare the temperature dependence of 
p s (T) for the aged sample and the BdG response cal- 
culation in the limit of strong impurity scattering. The 
impurity concentration nj mp = 6% was chosen to repro- 
duce the observed suppression of T c and p s {0) of the aged 
vs. fresh sample. The relevant parameters are given in 
the caption to Fig. [7J The BCS coherence length is given 
by £o ~ ftv.F/(7rAo), where the Fermi velocity is given 
approximately by vf ~ aot/Ti with Aq = 0.187t, so that 
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FIG. 7: Temperature dependence of the normalized super- 
fluid density p s (T)/p s (0) = \(0) 2 /\(T) 2 in PuCoGa 5 . The 
solid squares are for Ho = 60 mT and the open squares for 
H = 0.30 T after 400 days, using cr v (0) = 0.45(1) ps' 1 cor- 
responding to A(0) = 487(10) nm. The solid line is the short- 
coherence-length BdG superfluid density (for p s (0;ni mp = 
6%)/p s (0; 7ii mp = 0%) ~ 0.35) in the strong scattering limit 
described by Ao = 0.187t, T c — 0.115£, and impurity concen- 
tration nimp = 6%. For the pure case, ni mp = 0%, one finds 
Ao = 0.314£ with T c o = 0.143£. The statistical error bars are 
due to averaging over 20 random impurity configurations at 
fixed concentration. 

£o ~ 2ao- The excellent agreement between the short- 
coherence-length BdG calculation and the measured su- 
perconducting properties of the aged sample, and, at 
the same time, the failure of the dirty d-wave calcu- 
lation within AG pair-breaking theory to describe the 
data, indicates that the uniform, dilute-impurity theory 
by Abrikosov and Gor'kov is not applicable to PuCoGa5, 
where superconductivity is seemingly not uniformly sup- 
pressed. Specifically, both the BdG and AG dirty d-wave 
models require kfl 3> 1, where kf is Fermi wave vec- 
tor. However the former case requires strong impurity 
scattering with on-site suppression of Ao and fc/£ ~ 1, 
while the AG pair-breaking theory requires a uniformly 
suppressed Aq with £/Z < 1. 



V. DISCUSSION AND SUMMARY 

We now summarize our principal results. The radioac- 
tive decay of 239 Pu allows one to study the effects of pair- 
breaking defects on the properties of the superconductiv- 
ity in PuCoGa 5 . The low-temperature, quasi- linear tem- 
perature dependences of the superfluid density and pen- 
etration depth in both fresh (25 day-old) and aged (400 
day-old) PuCoGas are consistent with a line of nodes in 
a d-wave order parameter. 

The fresh sample is almost defect free, as evidenced by 
the small reduction in T c , e.g., AT C /T c0 w 1.4%. We 
have, therefore, compared our data in the fresh sam- 
ple to a dirty d-wave model for a strong-coupling su- 
perconductor in the presence of strong impurity scat- 
tering. We find that the zero-temperature penetration 
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depth and superconducting gap are A(0) = 265(5) nm 
and A = nkgT c0 with n — 3.0(1), respectively. The lat- 
ter is enhanced compared to the weak-coupling d-wave 
case, where n — 2.14. The model scattering rate used 
was consistent with the slight reduction in T c for the fresh 
sample, according to conventional AG pair-breaking the- 
ory. Our results on the fresh sample are consistent with 
recent NQR and NMR experiments in fresh PuCoGas 

In the case of the aged sample, we argued that the 
presence of strong pinning in PuCoGas allows an accu- 
rate measure of the temperature dependence of A(T), but 
that distortion of the FLL (which broadens the linewidth) 
means that the measured magnitude of A(0) is only a 
lower limit. The dirty d-wave model can also fit the tem- 
perature dependence of the superfluid density p s (T) (or 
A(T)) with A reduced by about 20%, but with a scat- 
tering rate nearly an order of magnitude smaller than 
predicted by conventional AG pair-breaking theory. Fur- 
thermore, the dirty d-wave model is unable to account 
for the reduction of at least 70% in p s (0). This is at- 
tributed to the fact that PuCoGas possesses a relatively 
short coherence length, and, therefore, the conventional 
AG pair-breaking theory, in which the order parameter 
is spatially averaged over impurity sites, is inappropri- 
ate. This was pointed out by Franz et ali^ for damaged 
or doped YBCO superconductors, but in that work a 
comparison with the data was made only for zero tem- 
perature. 

Accordingly, we have modeled the full temperature de- 
pendence of p s (T) in aged PuCoGas by solving the BdG 
equations for a short-coherence-length, weak-coupling d- 
wave superconductor in the strong scattering limit. (For 
small scattering rates the short-coherence-length model 
agrees fairly well with the dirty d-wave modeli^i) The 
BdG model is able to reproduce quite well the tempera- 
ture dependence and magnitude of p s for a nominal 6% 
impurity concentration, chosen to reproduce the magni- 
tude of p s (0). This impurity concentration is reasonable, 
though we attach no particular importance to its exact 



value. Like the dirty d-wave model, the BdG model is 
consistent with a reduction of Ao in the aged material, 
compared to the fresh sample. Both the dirty d-wave 
and BdG models are consistent with A(0) = 498(10) nm 
in aged PuCoGas, where we remind the reader that be- 
cause of FLL distortion this is a lower limit . 

Computational feasibility dictates that the BdG lat- 
tice model calculations are for a two-dimensional sys- 
tem. Therefore, the derived gap parameters are only 
semi-quantitative. (The values of A(0) are accurately de- 
termined from Eq. (J3J), however, and our quoted values 
are derived from this expression.) Nevertheless, we do 
not expect strong deviations between 2D and 3D mod- 
els for the in-plane penetration depth, and, therefore, 
believe that a 2D model should provide a good qualita- 
tive description of the essential physics in this interesting 
superconductor^ 

Our data and analysis in aged PuCoGas suggest that 
although small parts of the sample are significantly disor- 
dered, and consequently the superfluid density is strongly 
suppressed in these regions, superconductivity remains 
remarkably resilient. 
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